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§ . ABSTRACT 
On 

The MFA approach for simulations with dynamical fermions in lattice gauge 
Qh! theories allows in principle to explore the parameters space of the theory (e.g. 

the (5, m plane for the study of chiral condensate in QED) without the need 
of computing the fermionic determinant at each point. 

We exploit this possibility for extracting both the renormalization group tra- 
jectories ("constant physics lines") and the scaling function, and we test it in 
the Schwinger Model. We discuss the applicability of this method to realistic 
theories. 

In this paper we present a further test of the potentialities of the Microcanonical 
Fermionic Average (M.F.A.)|l| method. We apply this procedure to the Schwinger 
Model. Since in this approach, the main computer cost resides in the evaluation of an 
effective fermionic action at fixed pure gauge energy by evaluating all the eigenvalues 
of the fermionic matrix at m = 0, it is possible, at no extra cost, to move in the plane 
P, m. 

We have made full use of this fact for obtaining the constant physics trajectories 
of the model, whithout any a priori assumption of their form, which is exactly known 
in the Schwinger Model, since here the Renormalization Group amounts to simple 
dimensional analysis. 

The continuum Schwinger model || is confining; it is exactly solvable at zero 
fermionic mass, so that we can compare the results of our simulations with exact 
ones. It is superrenormalizable and the Callan-Symanzik (3 function vanishes. The 
chiral current is anomalous; if the chiral limit is obtained from m ^ 0, then the 



6 = vacuum is selected. In this vacuum the chiral condensate is (with one flavour) 
e~ 1 (ipil ! ) c = 0.15995 This is the value of the chiral condensate to be compared with the 
results of lattice simulations, where its chiral limit is obtained from m^O. 

The M.F.A. method is fully described in JT[|. Starting from the partition function 
written as sum of the fermionic and pure gauge contributions, we define the density of 
states at fixed pure gauge Action (i.e. Euclidean Energy) and an effective fermionic 
action which is (minus) the logarithm of the average of the fermionic determinant over 
fixed energy configurations. All the thermodynamical quantities can be computed as 
1-dimensional integrals; bosonic operators can be directly computed at m = In the 
present simulation we use 1 flavour of staggered fermions, and the pure gauge part is 
described in terms of non compact fields, so that the density of states is analytically 
known []. 

Since e c is dimensionful, 8 explicitely contains the lattice spacing: 8 = so that 
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the continuum limit of the theory is approached, assuming super-renormalizability, at 
8 — » oo. The limit must be reached keeping fixed the dimensionless ratio ^ = \f]3m. 
This ratio defines constant physics trajectories. We present here the results for the 
100 2 lattice. 

We compute all the eigenvalues of the fermionic matrix; this allows us to compute 
exactly the fermionic determinant, and hence the Effective Action for all values of the 
mass, including m = 0. We have done so for 8 values of the energy, from 0.02 to 0.25. 

The chiral condensate at m = vanishes on a finite lattice, so it must be obtained 
as the limit m — > 0. To reach the correct continuum value, this limit has to be taken 
simultaneously with the — > oo one, keeping the product y//3m fixed. This can be 
easily done with this method, which does not require a separate simulation of the 
fermionic contribution for each pair of parameters (8,m) [0]. 

We have repeated this procedure for 10 values of m c /e c , always finding a clear 
scaling window, where we compute the chiral condensate. Using a quadratic fit for the 
experimental points, we obtain for the intercept: e c T 1 (^^) = 0.159 ± 0.003 in perfect 
agreement with the theoretical value. 

To better show the potentialities of the method, we have investigated the possi- 
bility of defining constant physics trajectories without any a priori knowledge of their 
behaviour. This procedure would be the appropriate one in any perturbatively renor- 
malizable theory like QCD. 

We identify as physical observables the continuum chiral condensate and the 
charge, and look for constant physics trajectories in the plane (0,m). 

Since a{ipip)c = (i'^L, and that the plaquette behaves as 1/0 as 8 — > oo: in this 
limit (E)l — > e 2 c a 2 . This implies that 



is independent on the cutoff and then defines the constant physics trajectories as the 
values in the 8,m plane where R is constant. 
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Constant level lines on the surface defined by R are thus the trajectories we are 
looking for, and they are reported in Fig. 1, compared with = constant. 

It can be seen that the agreement is indeed very good, apart from large masses 
and/or small j3 where we are affected by finite cutoff effects. 

It is important to discuss also the applicability of these results to more realistic 
theories, for instance QCD. The choice of the operators to be fixed in the renormal- 
ization procedure is dictated by physical motivations, and the choice of the continuum 
charge is intrinsic of the Schwinger model. 

On the other hand, to fully exploit the potentialities of the MFA method, one 
is restricted to the use of either bosonic operators, or fermionic ones related to the 
eigenvalues of the fermionic matrix. 

In QCD a convenient choice would be the ratio of the longitudinal to transverse 
susceptibilities || R(f3,m) = Xl/xt where both numerator and denominator are com- 
putable from the eigenvalues of the fermionic matrix. 

Once constant physics trajectories (which are renormalization group trajectories 
in renormalizable theories) have been defined, any other physical observable should 
be constant if computed along them. For instance, if we compute hadronic masses in 
lattice units m^a, all their dependence on bare parameters along R.G. trajectories will 
come from the lattice spacing a. In this way one can extract the scaling functions. 
This allows to extract the continuum values of the observables. 

This work has been partly supported through a CICYT (Spain) - INFN (Italy) 
collaboration. 



FIGURE CAPTION 

Constant Physics trajectories in (/3,m) plane (continuous lines); circles are /3 1//2 m = 
constant lines. 
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